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The use of reverse time chaos allows the realization of hardware chaotic systems that can operate at speeds 
equivalent to existing state of the art while requiring significantly less complex circuitry. Matched filter 
decoding is possible for the reverse time system since it exhibits a closed form solution formed partially by 
a linear basis pulse. Coefficients have been calculated and are used to realize the matched filter digitally 
as a finite impulse response filter. Numerical simulations confirm that this correctly implements a matched 
filter that can be used for detection of the chaotic signal. In addition, the direct form of the filter has been 
implemented in hardware description language and demonstrates performance in agreement with numerical 
results. 
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Recent developments in the field of chaos have 
shown that certain chaotic systems exhibit closed 
form analytic solutions. Because of the unique 
form of these solutions, it has been shown that 
a matched filter can be designed to detect the 
presence of chaotic waveforms generated by these 
exactly solvable systems. Hardware has been de¬ 
veloped realizing both the chaotic system (in the 
form of a chaotic oscillator) and the matched fil¬ 
ter, but both of these have been limited to rel¬ 
atively low frequencies. Reverse time chaos has 
been presented as a means to increase the achiev¬ 
able operating frequency by redefining the chaotic 
system equations in a form that allows a hybrid 
digital and analog circuit realization. This work 
expands on reverse time chaos by deriving a nu¬ 
merical model of a matched filter for the reverse 
time system. In addition, a digital FIR imple¬ 
mentation of this filter is demonstrated to func¬ 
tion correctly in simulation. 


I. INTRODUCTION 

Although many types of chaos continue to evade ana¬ 
lytic solutions, the assumption that no solvable chaotic 
dynamics exist has been broken by the discovery of hy¬ 
brid (consisting of interdependent continuous time and 
discrete time states) chaotic systems that exhibit exact 
solutions^ 3 . A linear matched filter, which can be used 
to exactly recover the discrete-time state of the chaotic 
system when driven with only that system’s continuous¬ 
time state, has also been developed using these solu¬ 
tions. Matched filter decoding significantly reduces the 
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difficulty of integrating chaos into a wide range of elec¬ 
tronic systems, including both wireless communication 
and radar. 

One of the most promising areas of research in 
chaotic systems has long been in their utility for 
communication^. Previous work has shown that such 
systems can be constructed and tuned such that they 
possess advantageous dynamics for various communica¬ 
tion schemes^. In addition to possessing these tar¬ 
geted dynamics, chaotic systems can also be subjected to 
control schemes that can maintain desired trajectories^. 
By combining these characteristics with matched filter 
decoding^, many of the necessary components for a 
modern high-performance communication system may be 
realized with chaos-based systems. 

Various elements of such systems have been analyzed 
to determine how they might perform in a number of 
applications. In the presence of additive Gaussian white 
noise, chaotic systems can exhibit error rates comparable 
to communication systems currently in us e^S>U . Many 
issues associated with multipath propagation may also be 
mitigated when the propagated waveform is chaotioiS. 

Reverse time chaos provides a potential solution for 
realizing chaos in hardware with both solvable and con¬ 
trollable properties without sacrificing the ability to scale 
in speed. First proposed by Corron et al. in~, reverse 
time chaos describes behavior that differs from tradi¬ 
tional chaos by using the current state of the system 
to represent all of its past states instead of all of its fu¬ 
ture states. Despite this difference, reverse time chaos re¬ 
tains a positive Lyapunov exponent and a corresponding 
sensitivity to initial conditions that defines traditional 
chaotic systems. The many similarities between tradi¬ 
tional (or forward time) and reverse time chaos also al¬ 
low for formulation of an exact solution. Because of this, 
in a method analogous to that used for the forward time 
systems in^, the development of a linear matched filter 
that can be used to detect the original chaotic signal is 
also possible. 
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In this paper, we build on recent work demonstrating 
a functional reverse time chaotic oscillator™ by devel¬ 
oping an expression for its corresponding matched filter. 
Realization of the matched filter completes the transmit¬ 
ter and receiver system presented in previous work using 
forward time chaos discussed above. With this complete 
system, many of the advantages cited for the forward 
time case can also be realized with reverse time chaos. 



Figure 1. Numerical simulation of exactly solvable system in 
shift band (time evolution) 


II. EXACTLY SOLVABLE CHAOS 


A. Forward Time 


The chaotic system used by Corron et al. can be de¬ 
scribed by a set of two equations written in terms of two 
states. These consist of a continuous time state u(t) G R 
that provides exponential growth so that the system does 
not die out and a discrete time state s(t) € {—1,1} or 
s(t) G {0,1} that maintains boundedness so that the sys¬ 
tem does not grow unbounded. In order to differentiate 
its behavior from the reverse time system discussed in 
the next section, the dynamics described by this behav¬ 
ior are those that have been previously referred to as for¬ 
ward time chaos. u(t) is defined by the continuous-time 
differential equation: 


u — 2/3u + (w 2 + f3 2 ) * (u — s) = 0 (1) 


where to = 2 tt and 0 < /? < In 2. Various system dynam¬ 
ics can be realized by selecting specific definitions of s(£); 
for Lorenz-like dynamics, s(t) is defined as: 


u = 0 —> s(t) = sgn(u(t)) (2) 


where 


sgn{u) 


— 1 u < 0 
1 u > 0 


(3) 


When operating with this s(t), the system is said to be 
in the shift band. Figure [j] demonstrates this behavior 
graphically. 


As with many other chaotic systems, the phase por¬ 
trait of this system provides more insight into long term 
dynamics. The phase portrait for the shift band is pre¬ 
sented in Figure 0 



Figure 2. Numerical simulation of exactly solvable system in 
shift band (phase portrait). 


Although the system dynamics are captured by u and 
u, s is also included to create a three-dimensional plot; 
this allows for an alternate visualization of the system 
oscillations u being centered around the discrete levels of 
s. Despite the complexity this system appears to pos¬ 
sess when viewed graphically, it can be shown to have an 
analytic solution given as: 

OO 

u(t) = ^2 s mP{t - ™) (4) 


{ (1 — e /3 )e /3t (cos (cot) — ^ sin(wt)) t < 0 
1 — e /3 ( t_1 ^(cos (ojt) — ^ sin (uit)) t = 0 (5) 

0 t > 0 

where P(t) acts as a basis function and s m G —1,1 is 
the symbol sequence defined for all Although the 
simple existence of this solution is remarkable, its most 
significant characteristic is in its description in Equation 
|4] in terms of convolution of a basis function with a ran¬ 
dom bipolar pulse. This description notably allows for 
derivation of the matched filter that can be used to par¬ 
tially regenerate s(t) when presented with only u(t). The 
matched filter can be described by: 

e + 2^ + (w 2 +/3 2 )£ = (o , 2 +/3 2 )v(t) (6) 
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where 

f) = v(t + 1) — v(t) (7) 

In Equations [6] and [3 £ represents the output of the 
matched filter (which is also the recovered s(t)) while 
v(t) represents the input, which is assumed to be u{t)^. 
Previous work has demonstrated that hardware imple¬ 
mentations of this filter can be used to successfully re¬ 
cover the original s(i)I&. 

B. Reverse Time 

Although both components of the forward time exactly 
solvable system have been extended to circuit simulations 
operating near 2 published hardware implemen¬ 

tations operate in the audio frequency rang o 2 ’ 18 . Each 
circuit consists of many analog components that scale 
poorly with frequency and do not integrate well with 
modern high density integrated circuit (IC) processes. 
In particular, the use of a negative impedance converter 
in many of the published circuits provides a significant 
challenge for high frequency designs of the oscillator and 
the requirement for a long delay time adds significant 
complexity to the matched filter. 

Reverse time chaos provides a method for generating 
functionally equivalent waveforms with significantly sim¬ 
pler circuit implementations than have been published for 
the forward time case2,ilZ. In contrast to these circuits, 
which use many active components that may be difficult 
to scale in frequency, the reverse time chaotic oscillator 
requires only three passive components (forming a se¬ 
ries resistor, inductor, and capacitor - RLC - tank) in its 
simplest form; the remainder of the circuit can be imple¬ 
mented digitally^. Complementing this, the reverse time 
matched filter presented in this work can be implemented 
entirely in digital circuitry. 

A mathematical description of reverse time chaos can 
be defined by the second order system: 

u + 2f3u + (w 2 + /3 2 )u = ( ui 2 + j3 2 )s(t) (8) 

as first described ini^. Both /3 and oj are control param¬ 
eters for this system; however, to maintain consistency 
with typical engineering practice, ui will be defined as the 
radian frequency of operation. s(t), the forcing function, 
can be described as 

s(t) = As n (9) 

where s n E {—1,1} is a random sequence with amplitude 
A , which can take any real value. 

Upon inspection, it is apparent that Equation [5] de¬ 
scribes a harmonic oscillator with positive damping pro¬ 
ducing oscillations centered around the value of s(t). An 
illustration of this behavior generated with numerical 
simulation for /3 = In(2) and ui = 2n is shown in Fig¬ 
ure [3] 



Figure 3. Numerical simulation of reverse time system in shift 
band (time evolution). 


The positive damping term causes the value of u(t) 
to decrease in magnitude and is the reason that the re¬ 
verse time system requires s(t) to act as a forcing func¬ 
tion. Without this forcing, u(t) would simply die out 
after ringing down to 0. 

As described in the introduction, a reverse time chaotic 
system differs from a forward time chaotic system by us¬ 
ing its initial condition to describe the set containing all 
of its previous states rather than its future states. To 
better visualize this, Figure [3 below shows Figure |3] with 
its x axis reversed. 
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Figure 4. Numerical simulation of reverse time system in shift 
band (time reversal). 


When viewed in this manner, the dynamics of u(t) ap¬ 
pear identical to those of |T] operating in the shift band 
(Figure [[]). Further confirmation of the similarity is 
shown by examining the phase space projection, shown 
in Figure [3 
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Figure 5. Numerical simulation of reverse time system in shift 
band (phase portrait). 

The reverse time system’s phase portrait matches al¬ 
most exactly with the phase portrait from the forward 
time system apart from the levels of s around which the 
trajectories orbit. In the previous section, these were - 
1 and 1, while in Figure [5] they are -2.5 and 2.5, which 
corresponds to the value of A used in the reverse time nu¬ 
merical analysis routine. As in the forward time system, 
a solution can be described in terms of a basis pulse; the 
basis pulse for the reverse time system is given as: 


Ug(t) 


f° 


1 l-e-' 34 
oj' z +(3' 2 

(cos(wf) + (^) sin(wf)) 

e-* 

(cos(wf) + (^) sin (ut)) 



t < 0 
0 < t < 1 
t > 1 


( 10 ) 

A plot of Equation [TUI evaluated for 0 < t < 3 is shown 
in Figure OH 
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Figure 6. Reverse time basis pulse. 

The difference between the two cases in Equation ITOl is 
clearly visible; for 0 < t < 1, the basis pulse is offset by 1 
while for t > 1, this offset is removed. As expected, this 
response is identical but reversed in time to the forward 
time basis pulse in Equation [5j Because the system’s 
solution can be described in terms of this basis pulse, it 
can be extended to all 

t 

u(t) = ^2 s n Ug(t-n ) (11) 

n =—oo 

The existence of this solution also allows for the devel¬ 
opment of a matched filter with the same process used 
for the forward time system. 


III. MATCHED FILTER 


A common problem in many signal processing appli¬ 
cations is to determine if a known signal is present in a 
received signal that may contain significant noise. Al¬ 
though there are multiple approaches to solve this prob¬ 
lem, the use of a filter whose response is matched to the 
known signal - a matched filter - is frequently the so¬ 
lution of choice. This matched filter can be described 
as the system with the maximal output SNR when pre¬ 
sented with a known signal corrupted by additive white 
Gaussian noise (AWGN). A brief mathematical descrip¬ 
tion fellow s hips. 


A. Derivation 


Given a known transmitted signal x(t), the received 
signal v(t) can be written as: 

v(t) = Ax(t — t p ) + n(t) (12) 


where A is amplitude scaling due to attenuation or gain 
of the channel, t p is the propagation delay through the 
channel, and n(t) is AWGN added in the channel with 
a power spectral density (PSD) S n (u>). The SNR of this 
system can be defined as: 


SNR = 


max (| Ax(t — to)| 2 ) 
|n(t) 2 | 


\Ax(0)\ 2 

\n(t )\ 2 


(13) 


The matched filter is the linear filter {h(t) or H(co)) 
that maximizes its output signal-to-noise ratio (SNR) at 
some instant relative to a starting point to and can be 
can be written ash; 


h(t) = kx(t x — t ) (14) 


Equation [14] describes a filter whose response is identi¬ 
cal to the known transmitted signal x(t) but reversed in 
time and shifted by some amount t x (k may be used to 
adjust the amplitude, but typically is left equal to 1). For 
the reverse time system, the received response is taken to 
be the basis pulse given by Equation [TO] The resulting 
required matched filter response is shown in Figure [7] for 
/? = ln(2). 
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Figure 7. Reverse time matched filter response. 


B. Verification 

Using a numerical model for the matched filter, its per¬ 
formance in a number of scenarios of interest was eval¬ 
uated. The first test performed was to determine if the 
filter could properly detect the basis pulse when inserted 
into a random signal. Figure[5]demonstrates this process. 
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Figure 8. Matched filter response to basis pulse, (top) Input 
to filter and (bottom) output from filter. 


The input (v(t)) is shown in the top while the output 
from the numerical filter (v 0 (t)) is shown in the bottom 
plot. These results indicate that it correctly responds to 
the presence of a basis pulse. To verify the expected re¬ 
silience in the presence of AWGN, three additional tests 
were performed with increasing noise levels added for 
each successive run. MATLAB’s awgn function was used 
to add noise to the existing model. Figure [9] shows these 
results with the SNR set to 0 dB. In this figure, the seg¬ 
ment containing the basis pulse with added noise is high¬ 
lighted in red. 
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Figure 9. Matched filter response to basis pulse in presence 
of AWGN (SNR=0 dB). (top) Input to filter with basis pulse 
overlaid in red and (bottom) output from filter. 

Despite the significant noise added to the first instance 
of the basis pulse, the matched filter was able to correctly 
detect its presence with approximately the same response 
magnitude as generated by the basis pulse alone. The re¬ 
sponse to the random segments at the beginning and end 
of the input signal remained low throughout; this result 
combined with the peak behavior previously discussed 
indicates a strong selectivity for the basis pulse. 

C. Chaotic Input/Matched Filter Output 

While the ability to detect a single basis pulse pro¬ 
vides verification that the matched filter model is oper¬ 
ating correctly, its true utility comes from the ability to 
detect a reverse time chaotic signal. This signal consists 
of many instances of the basis pulse summed together, 
therefore the matched filter output shows a large num¬ 
ber of peaks as the input is shifted through its response. 
The magnitude of these peaks is diminished as compared 
to the response when presented with a single instance of 
the basis pulse, but will be shown in the next section 
to be high enough to allow for reliable detection of the 
chaotic waveform. Reusing the numerical model from 
the previous section, a chaotic signal was inserted into a 
random waveform to verify the filter’s functionality. The 
output of this test is shown in Figure ITOl 



t 

Figure 10. Matched filter response to reverse time chaotic 
waveform, (top) Input to filter and (bottom) output from 
filter. 

For completeness, the effect of AWGN on the chaotic 
waveform has also been explored with the numerical 
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matched filter. The same noise level (SNR = 0 dB) used 
peviously was selected for this analysis. These results are 
presented in Figure fill 



Figure 11. Matched filter response to reverse time chaotic 
waveform in presence of AWGN (SNR=0 dB). (top) Input to 
filter with chaotic signal overlaid in red and (bottom) output 
from filter. 
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Figure 12. Threshold levels used for reconstruction of s n - 
(top) Input to filter and (bottom) output from filter with 
threshold levels indicated. 

New values of s n always correspond to the matched 
filter output going above (or below) the high and low 
threshold values. Simply detecting these events, how¬ 
ever, does not properly reconstruct s n - it only indicates 
the presence of the basis pulse (it is important to re¬ 
iterate that the matched filter is matched to the basis 
pulse, not the overall chaotic waveform or s n ). The so¬ 
lution in Equation |TT] allows for s n to hold its value for 
many cycles, which results in many oscillations around 
the threshold. To account for this, the midpoint thresh¬ 
old is used to detect when a switch in the s n value has 
occurred and then update the recovered s n . This process 
is shown graphically in Figure [T51 


As with the basis pulse, even high magnitude AWGN 
does not significantly diminish the response of the 
matched filter to the presence of a chaotic waveform. 


co 
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D. s„ Reconstruction 


Limited post-processing can be applied to the matched 
filter response to reconstruct the s n sequence used to gen¬ 
erate the chaotic signal applied to its input. This post¬ 
processing method requires the establishment of three 
threshold levels: a high threshold, a zero crossing or mid¬ 
point threshold, and a low threshold. These thresholds 
were determined by visual inspection in this work; al¬ 
though non-optimal, this method still allows for proper 
reconstruction. Figure [12] demonstrates the placement of 
these levels where the high and low thresholds are de¬ 
noted by the upper and lower red lines and the midpoint 
threshold is denoted by the black line. 


Figure 13. Reconstructed s„. (top) Original s signal from os¬ 
cillator and (bottom) reconstructed s from filter output with 
each decision point indicated by a blue asterisk. 

When compared to the original s state of the reverse 
time chaotic oscillator (top), the reconstructed s (bot¬ 
tom) appears identical apart from the time shift intro¬ 
duced by the matched filter operation. The blue asterisks 
marked on the reconstructed s plot indicate the decision 
made by the post processing algorithm at each time step. 
For the vast majority of the steps, this value is 0, corre¬ 
sponding to no change in s (i.e., many asterisks are re¬ 
peated at reconstructed s = 0). When the matched filter 
output crosses through the midpoint threshold and then 
exceeds either the high or low threshold, however, the 
decision changes to indicate a new value. These decision 
points are then used to generate the overall reconstructed 
s as shown. 
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IV. FIR IMPLEMENTATION 


Numerical modeling in MATLAB provided a straight¬ 
forward path for realization of the matched filter in hard¬ 
ware as a digital FIR filter. The selection of a digital fil¬ 
tering method provides two distinct advantages: it allows 
for a highly precise definition of the filter response while 
also side stepping the challenges raised by attempting to 
implement the inherently unstable reverse time matched 
filter response in analog hardware. Values used for the 
filter’s magnitude response can be used directly to define 
the coefficients necessary for the FIR algorithm. Any 
N th order constant coefficient digital FIR function can 
be described as a convolution sum as shown in Equation 

H3-1. 


x[n] 


I 



Figure 14. FIR filter block diagram. 


N-l 

Y[n] = X[n] * H[n] = H[k] * X[n — fc] (15) 

k=0 


The exact function performed is determined by the val¬ 
ues of H, which are the FIR coefficients. This form can 
be implemented in hardware directly (and is often re¬ 
ferred to as the direct-form FIR), but is generally less 
desirable due to its use of multipliers. One solution for 
addressing this shortcoming is using the sum of power of 
two (SOPOT) method to decompose multiplication into 
addition and multiplication by powers of two 22 ’ 23 . Equa¬ 
tion [TS] describes the process of determining the SOPOT 
representation of an arbitrary value s(n). 

J -1 

s(n) = ^2 a i,n * 2 bi,n (16) 

i =0 

where ai tTl G {—1,0,1}, G (0,1,..., u}, and J is the 
number of terms necessary to represent s(n )22. The max¬ 
imum value that can be represented by SOPOT is related 
to u, so its value determines the allowable range of co¬ 
efficients. Because SOPOT contains only multiplication 
by powers of two, the multiplication functions can be re¬ 
place by logical left shifts that shift by the bi. n values. An 
algorithm based on the work in^i was developed in MAT- 
LAB to automate the calculation of these coefficients for 
any real number. 


A. Architecture 

A general block diagram of the FIR architecture is 
shown in Figure [TTP 5 . 


From this figure, it is apparent that the algorithm re¬ 
quires relatively few components for its hardware real¬ 
ization. The unit delay blocks (indicated as z _1 ) map 
directly to D flip-flops where the clock period sets the 
delay time. The weights (denoted by /[0] to f[N — 1]) 
require a varying number of adders and left shift circuits 
as indicated by the SOPOT decomposition. A MATLAB 
routine was also written to generate and combine these 
blocks automatically, given a set of FIR weights. 


B. HDL Simulation 


Making use of the Verilog code generated for the FIR 
filter, a HDL simulation was performed using ModelSim 
to determine expected performance of a physical imple¬ 
mentation. 32 bit values were used for internal calcula¬ 
tions to allow for high precision while 10 bit values were 
used for the input to correspond to typical ADC resolu¬ 
tions. An operating clock frequency was set to 180 MHz 
due to the use of the 1.8 MHz chaotic oscillator with 100 
FIR coefficients. A simulation, shown in Figure HU con¬ 
firms the code’s performance in the presence of AWGN. 
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Figure 16. Simulated BER for basis pulses of varying length. 

From the above results, it is immediately apparent that 
the best performance is realized with t = 3. This result 
was selected for a comparison against the previously pub¬ 
lished theoretical performance^. Results of the compari¬ 
son are shown in Figure [TT] 


Figure 15. Simulation results of FIR matched filter with 
AWGN (SNR=0 dB). (top) Input to filter with chaotic signal 
overlaid in red and original s signal, (middle) output from 
filter with threshold levels indicated, and (bottom) recon¬ 
structed s from filter output with each decision point indi¬ 
cated by a blue asterisk. 


As expected, the Verilog model is able to correctly de¬ 
tect the chaotic signal despite the significant magnitude 
of added noise. The large scales used for the y axis on 
both the filter input and output result from converting 
the binary values generated from the simulation into dec¬ 
imal directly. 


C. System Performance 



Figure 17. Comparison of theoretical and simulated BER. 

While the theoretical BER takes in to account the error 
induced by the determinism inherent in a chaotic wave¬ 
form, it does not address addition errors introduced by 
the digitization and truncation of the matched filter that 
was performed in this work. As expected, with these in¬ 
cluded, the performance of the overall system is slightly 
reduced as compared to theoretical expectations. 


With the reconstructed s sequence, performance of 
the oscillator and matched filter combination can be ex¬ 
plored. A common method of evaluating communication 
systems is to examine their bit error rate (BER), which 
can be described by: 

number of incorrect bits , 

BER =-^---—- (17) 

number of total bits 


Evaluating the BER. for varying lengths (maximum 
values of t) of the basis pulse in Equation [TO] allows for 
selection of the optimal length. Figure [E)] demonstrates 
the effects of varying t from 3 to 5. 


V. CONCLUSIONS 

Matched filter decoding has been shown to provide a 
viable means of detection of reverse time chaotic wave¬ 
forms. Building on previous work, a numerical model 
was constructed and analysis predicted that these wave¬ 
forms could be reliably detected even in the presence of 
high magnitude AWGN. Both this ability and the abil¬ 
ity to reconstruct the original s sequence using a simple 
post-processing scheme have been demonstrated. 

Conversion of the numerical model to a FIR filter has 
been performed and implemented in direct form. To more 
readily translate to a synthesizable design, this form was 
decomposed using the SOPOT method and translated 
into HDL code. Simulation results of this code indicate 
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that a hardware realization is straightforward and prac¬ 
tical on both a low cost FPGA and a high performance 
ASIC. 

Finally, these results confirm that reverse time chaos 
can be used in the same transmitter/receiver model as 
forward time chaos without requiring the complex and 
difficult to scale hardware associated with existing for¬ 
ward time chaotic systems. Applications in many areas 
including communications and radar can be satisfied with 
the simple hardware required for reverse time chaos. 
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